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Abstract 

The coherence between quantum states with different particle numbers — the Fock-space coher- 
ence — qualitatively differs from the more common Hilbert-space coherence between states with 
equal particle numbers. For a quantum dot with multiple channels available for transport, we find 
the conditions for decoupling the dynamics of the Fock-space coherence from both the Hilbert- 
space coherence as well as the population dynamics. We further find specific energy and coupling 
regimes where a long-lived resonance in the Fock-space coherence of the system is realized, even 
where no resonances are found either in the populations or Hilbert-space coherence. Numerical 
calculations show this resonance remains robust in the presence of both boson-mediated relaxation 
and transport through the quantum dot. 

PACS numbers: 72.10.Bg, 73.21.La, 73.23.-b, 73.23.Hk 
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I. INTRODUCTION 



According to quantum theory, a physical system can exist in a coherent superposition 
of distinct values of a single physical observable. One such superposition is over the num- 
ber of particles present in the system. Such Fock-space coherence (FSC) has been recently 
and controllably demonstrated in bosonic systemsP^ The lifetime of these novel states may 
be considerably enhanced due to the near quenching of single-particle relaxation channels. 
This mitigation of both decay and decoherence is a prime requisite, for example, in the 
demonstration of scalable quantum information processing architectures. In Fermionic sys- 
tems such as charge carriers in a quantum dot, the influence of Fock states on the dynamics 
and transport properties has been investigated^ but not the coherence properties of the 
states themselves. Here we show that coherent fermionic Fock states in multilevel quan- 
tum dots can exhibit decoherence times at least an order of magnitude greater^ than the 
more conventional Hilbert-space coherence (HSC) between states with fixed particle num- 
ber. This estimate remains valid even in the presence of a strong perturbation such as 
resonant tunnelling transport through the system as well as coupling to a bosonic bath. 
Our fully non-Markovian calculations reveal a memory-kernel whose elements factorise into 
non-interacting blocks, one of which contains all of the Fock-coherence. Our results are 
independent of the number of tunnelling pathways available through the system, although 
more than one seems necessary, indicating that an interference effect may be responsible for 
the robustness of these states. Our results further demonstrate how novel quantum states 
can exhibit remarkably useful properties given careful — though not prohibitive — tuning of 
bias, level spacing, and barrier anisotropy of the system. 



II. MODEL 

Our model is that of a single quantum dot weakly coupled to biased source and drain 
semi-infinite leads, such that only a small number of electronic channels are available for 
transport. We also consider particle-preserving relaxation through a coupling to a bosonic 
reservoir (phonons, for example). The total Hamiltonian is given by 

H = -f^ leads + -^qd + -^boson + %T + (1) 
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where H\ eads = ± — eVfo)cf^^ <i s (ci) describes the source and drain leads as non- 

s(d) 

interacting fermion systems shifted by the bias voltage V\,. The creation and annihilation 
operators for the source (drain) leads are dJ, d y d s ^, respectively, and e s ^) are the respec- 
tive single-particle energies. The quantum dot in the non-interacting regime is given by 
H qd = '^^(hoji + eV g )c\ci, where c| and c, are system creation and annihilation opera- 

i 

tors, and the single-particle energies hwi are all shifted by the applied gate voltage, V g . 
For the sake of clarity we neglected Coulomb interactions throughout the paper, as these 
are not expected to qualitatively change the results beyond an energy-level renormaliza- 
tionPS The coupling to the electronic reservoirs is described by a tunnelling Hamiltonian, 
Ht = (T[dlci + h. c.) , where T-^ is an energy-independent tunnelling coefficient 

i,r=(s,d) 

for a particle tunnelling from the single-particle state \i) in the dot to the source (drain) 
reservoir^. We consider the qualitative effects of a boson-mediated relaxation on the FSC, 
and introduce a bosonic reservoir and its interaction with the quantum dot, described re- 
spectively by H hoson = y^ y eib\b h and H n = ^ L^cjc^ + h.c.j, where A ijt is a generic 

l i,j,l 

coupling coefficient containing the details of a specific electron-boson interaction. 



III. TIME EVOLUTION AND THE MEMORY KERNEL 

The non-Markovian time evolution of the quantum dot is considered in the weak coupling 
approximation (Born approximation), where terms up to second order in the interaction 
Hamiltonians Hj- and T-i^i are kept. 

The generalized master equation for the reduced density matrix^ is given by 

Pab(t) = V fdt' Pcd (t')r abcd (t - ay^\ (2) 

c,d J ° 

where the total memory kernel T a b c d(t) describing the dynamics of the system can be ex- 
plicitly derived^ directly from the microscopic Hamiltonian Eq. ([!]). Under a weak cou- 
pling approximation T a b c d(t) can be written as the sum of a transport-dependent (T) and 
relaxation-dependent (TZ) transition tensors, 

^abcd(t) = Tabcdif) + T^-abcdif)- (3) 
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The absence of bilinear terms in the interaction Hamiltonians Hf and Tin is due to the 
trace over the electronic reservoirs, where terms of the form (r\di\r) = (r\dj\r) = 0. 
The transition tensor T due to electronic transport is given by, 



T abcd (t) = £ {Ke [^sA^L ~ ^ T y(t)A a J ab _ 

%V (tf^abcd ~ n^V (*) A3L] } , (4a) 



with, 



fi£ y (t) = ( e ^«* - e^«*) /i, (4b) 
K^iN r (T r jTyh, (4c) 
Al d = (a|i |c)(d| C/! |6> - (a|c)(6|c Q 4|rf), (4d) 

and where the Latin indexes denote many-body system states, Greek indexes denote single- 
particle states, and r = s, d denotes the source and drain leads. Both leads are assumed 
to be metallic with occupied-state energies lying between a lower bound 4>b and the lead's 
chemical potential /j, s ^ d \ and unoccupied-state energies lying between the lead's chemical 
potential and an upper bound 4>t- The range between the bounds 4>b and 4>t are taken to 
be large relative to the characteristic energies of the system. Finally, N r denotes the density 
of states in lead r, and the frequencies = u>i — ujj. 

Similarly, the transition tensor due to boson-mediated relaxation is given by 

n abcd (t) = gJ2 i M ™* (eSe" ytt + ejSge-**'**) 

ijpqk 

+ M* jpqk (6gje*-«** + 8jge-*-«**)] (5a) 

where, 

= (a\c\c,\c)(d\c q cl\b) - {b\d){a\c\c pCj c^c), (5b) 
Mijpqk = AijkA pqk , (5c) 

and G is the density of states of the boson reservoir. 

IV. TWO TRANSPORT CHANNELS 

The transport and relaxation transition tensors presented in Eqs. (|4|[5]), are for a general 
number of available transport channels. For k transport channels, a minimum of 2 k dy- 
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FIG. 1. (a) Basis states for the two-channel case, defined as |0) = \N — l) gs , |1) = | A r ) gs , 
|2) = \N + l) gs , and |3) = |iV) es - (b) Representation of the non-zero memory kernel elements for 
the case of two transport channels. In the figure, the matrix elements represent the contributions 
to the dynamics between the reduced density matrix elements shown in the outermost row and 
column, due to transport (T), bosonic relaxation (1Z), or both. Not all elements of T a b c d are 
depicted; missing elements are the complex conjugates of those present. 



namical states are required, yielding 2 density matrix elements and 2 transition tensor 
components determining the dynamics of the system; despite symmetry-induced reductions 
in the number of independent components, the computational cost nevertheless increases 
exponentially with the number of channels. In what follows, we focus on the simplest non- 
trivial case of two available transport channels. 

The two-channel model is schematically depicted in Fig. fra. States |0), |1), and |2) 
respectively denote ground states with N — 1, N, and N +1 confined particles and state |3) 
denotes an iV-particle excited state. This is an experimentally accessible regime.^ The non- 
vanishing components of the memory kernel are shown in Fig. [T]d. The cell in the row labelled 
p a b and column p cd indicates whether the element T abcd is zero (empty cell) or depends 
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on a combination of transport (T) or boson- mediated relaxation (TV). For example, only 
transport contributes to T 3300 , only bosonic relaxation contributes to T 1133 , both transport 
and bosonic relaxation contribute to T 1113 , and T0022 vanishes. Elements not depicted in 
the figure may be inferred through the relation = Tl adc for a 7^ b and c 7^ d^ As 

expected, the only transitions where relaxation alone (TV) carries the dynamics are between 
states with the same number of particles and |3) in this case). 

In this matrix representation of the tensor T a b c d, two decoupled sub-matrices are ap- 
parent. The upper sub-matrix represents dynamics among population probabilities and 
coherence in states with identical particle number (HSC). The lower block describes coher- 
ence between states differing by one particle. Thus, the lower block describes the evolution 
of the FSC and is decoupled from the HSC. 

The decoupling of the FSC dynamics is not a result of the number of channels involved. 
For example, increasing the bias to allow more transport channels will increase the number 
of relevant elements of Tabd- The structure that emerges is similar to that depicted in 
Fig. [T)d, but the blocks themselves increase in both size and number; the FSC blocks remain 
decoupled from the HSC. Fundamentally, the structure is due to the fact that fermions 
cannot be individually created or destroyed through a scattering event (as opposed to bosonic 
systems). Since a bosonic bath does not alter the structure of the transition tensor, and 
since the time scale set by resonant tunnelling is generally much faster than that set by the 
bosonic coupling, we shall restrict ourselves to the transport problem and neglect the bosons 
in the remainder. 

A solution to the evolution equations, Eq. is obtained^ by moving to Laplace space 
and using the convolution theorem to transform the integrodifferential equations to a coupled 
set of algebraic equations. We numerically uncouple these equations and then transform back 
to the time domain by means of a Bromwich integral. 

The present transport model is described by seven parameters: Fermi energy Ep, bias 
voltage Vb, gate voltage V g , energy-level spacing AE = E a — Ep, orbital anisotropy e = 
1 — Tg/T* = 1 — Tfl/T^™ barrier asymmetry A = T^/T^ = Tg/TL and overall tunnel- 
coupling strength set by T^. The orbital anisotropy e describes differences in the tunnel 
coupling strength to different orbitals in the system. We define edge and core orbitals by 
the relative strength of their coupling to the reservoirs owing to the detailed shape of the 
wave function and the relative penetration into the tunnelling barriers.^ In relation to the 
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Parameter 


Ep 


eV B 


eV 9 


AE e 


A T* 


Value 


30 meV 


6 meV 


meV 


2 meV 1 


1 0.5 meV 



TABLE I. Standard set of parameter values used throughout the present work. Unless otherwise 
noted, these values were used for all results presented. 

2-channel case illustrated in Fig. [T^i, we denote orbitals |a) and \(3) as the edge and core 
orbitals respectively. 

In what follows, we plot our results as a function of one or two parameters while fixing 
the others to a set of "standard values." The standard values of the parameters are given 
in Table [TJ, which exhaustively describes our energetic configuration. For example, the two 
orbital levels have an energy of E a ^ = E F ± 4n + eV g . 

V. FSC RESONANCE 

We focus on the density-matrix element poi(£) which describes the FSC between the A^ — 1 
and the A^-particle ground states. (See Fig. [T^l.) The initial conditions for the density-matrix 
have an effect only on the amplitude of oscillations of the FSC, and not on the presence 
or location of a resonance. Thus, for the results presented here we choose a homogeneous 
initial FSC distribution as poi(0) = Po3(0) = Pi2(0) = P32(0) = 0.1. 

The time evolution of the FSC element poi(£) is shown in Fig. [2^, as a function of orbital 
anisotropy e (horizontal axis), where Tp = \T a with < % < 1. For all but a narrow range 
around e rcs ~ 0.22, the coherence drops off rapidly, within approximately 3 ps, corresponding 
to only a few tunnelling events (H/T a ~ 0.7 ps). This mirrors the decay found in the HSC,^ 
described, for example, by the matrix element pi3(t). Remarkably, at an orbital anisotropy 
of e Tes ~ 0.22, we observe a resonance in the FSC. The coherent oscillations persist for a far 
greater length of time. This resonance is unique to the FSC elements of p(t); they do not 
occur in either the HSC or the population probabilities p nn (t). To our knowledge, this is the 
first demonstration of a resonance in the Fock-Space coherence of an open quantum system 
of fermions. 

Figure [2]d shows the Fourier spectrum of the time-domain results shown in Fig. [2^,. At a 
given e, two frequencies dominate the spectrum. One is a high-frequency component which 
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FIG. 2. (a) Time evolution of the FSC element poi(t) as a function of the orbital anisotropy e. 
A resonance in the FSC is observed at an anisotropy of e res ~ 0.22 that maintains a persistent 
coherence well beyond the expected decoherence time, (b) Fourier transform of the results in (a). 
A linear dependence between frequency and e is evident. A spike in the low frequency amplitude 
occurs at the resonant anisotropy of e = 0.22. From this peak, we obtain the resonant frequency 
as given by Eq. ([7]). 

decreases linearly with e. This is the frequency related to the fast decay of the HSC, the 
population probabilities, and the FSC at very short times. There is also a low-frequency 
component evident in Fig. |2]d, corresponding to slow envelope oscillations. The resonance, 
clearly seen as the spike in Fig. [2]d, occurs on this low-frequency line. 

The resonant orbital anisotropy e res depends on the opacity of the source and drain 
barriers as well as on the overall strength of the tunnel couplings. Figure [3^, shows poi 
in the long-time limit (defined as the time beyond which the transients have decayed — 
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approximately 10 ps in the present function of both orbital anisotropy e and 

barrier asymmetry A. The evident peak (resonance) in the amplitude of poi a t this long- 
time limit gives the position of the resonant orbital anisotropy as a function of the barrier 
asymmetry. We observe a strong relationship between A and e rcs whose form is obtained 
from an empirical analysis (discussed later in this paper), and shown superimposed on the 
calculated data on Fig |3^l. A similar dependence is also obtained between the resonant 
orbital anisotropy and the overall tunnelling strength. In contrast, the dependence of e ms 
on the energetics is rather weak. Figure [3]d, for example, shows the FSC resonance, again in 
the long-time limit, as a function of orbital anisotropy e and level spacing AE = E a — -E73, 
and exhibits a linear and much weaker dependence of e res on AE. (Similar results are seen 
in the dependence of e res on the bias voltage Vb, and gate voltage V g .) 

As might be expected, the energetics of the system do play a role in the frequency of 
the coherent oscillations. Figure [5] plots the time evolution of the FSC element poi{t) as 
a function of the orbital level spacing AE for a resonant anisotropy e TCS = 0.22. The 
patterns seen in Figure [4] are suggestive of a coherent interference effect between tunnelling 
pathways as the source of the resonance. The strongest constructive interference for this set 
of parameters occurs for AE « 1.5 meV. 

We have empirically deduced that the dominant frequencies v evident in Fig. [2]d depend 
on the energy and coupling system parameters in the following way: 



i/±(e) 



A F 1 

E F + eV 9 — ± 21n(0 T )ro 2 (l + A 2 ) (1 - e) 



(6) 



We have further empirically deduced that the resonant frequency u Tes , also evident in Fig. [2}d, 
depends on the energetics of the system in the following way: 

v ies = ^(eV h + AE) + eV g . (7) 

At resonance, Eqs. (|6| and ([7]) describe the same frequency. Equating them yields the 
resonance anisotropy e res as a function of the Fermi energy, level-spacing, bias voltage, gate 
voltage, coupling strength, and barrier anisotropy: 

e -1 ( E F -AE-±eV h \ 

The fits to the numerical results presented in Figs. [3^, and (3jo were produced according 
to Eq. pi). Note that there are no free fitting parameters in this result. These result 
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FIG. 3. (a) FSC element poi in the long-time limit as a function of orbital anisotropy e, and 
barrier asymmetry A. The ridge denotes the resonance, (b) FSC element poi at the long-time limit 
as a function of orbital anisotropy and energy level spacing. The fits displayed in both plots where 
produced according to Eq. Q. 

may be used as an experimental guide when searching for this FSC resonance, as novel 
methods to extract electronic THz oscillations (equivalent to picosecond-range periods) are 
becoming accessible either by indirect electrooptical methods 1 18 * 19 ^ or by means of high-Tc 
superconductor devices^. 

VI. SOURCE OF FSC RESONANCE 

In order to address the source of the FSC resonance in our model, we consider the 
formation of a bound state between the quantum dot and the reservoirs EMU] g ucn a bound 
state may arise when the presence of a coupling term between the system and environment 
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FIG. 4. Time evolution of the FSC element poi(t) as a function of the level spacing within the 
bias window, and for a resonant orbital anisotropy e = 0.22. 

reduces the total energy. 

Denning the bound state of the QD and the reservoirs as = \QD)\R), where \QD) = 
P n \n) denotes the state of the QD, and \R) = ^ r f(ek)\ek)dek denotes the continuum 
of reservoir states for r = {s, d}, we solve the eigenvalue equation H\ip) = E\ip) for the total 
energy of the E, with the total Hamiltonian H given by Eq. [T} 

For a 2-channel case, the solution to the eigenvalue equation has the form,^, 

J2[E n 5 n>n ,- V / n ' [ n L de k ]P n , = EP n . (9) 

n'=0 r={s4} J ^ €k ~ E 

We note that the integration limits over the reservoirs become truncated by the single 
creation and annihilation operators of the tunneling Hamiltonian to J^, Fermi for removing a 
particle from, and ft, -T for adding a particle to, reservoir r. 

Fermi 

Therefore, the solution is found to be of the form, 

v /p\ TP T T \I<*AEi + -W,E 3 Pa,/3,E 3 + p / 1n \ 

Y {E) = E - l a , a ,E 3 ~ IhhEi. p p FT —J 1 = E, (lUj 

£j2 — £j — [i ai a,Ei + lp,p,E 3 \ 

with, 

Ia,b,E n — E n — E — ^ 
r={s,d} 

The existence of a bound state requires that Y(E) has at minimum one real solution for 
E < 0, which is equivalent to Y(0) < being satisfied. For the case of the transcendental 
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T r {T r 



Et, . e k 

t errm 



E 



de k + 



e k -E 



de k 



(11) 



equation, Eq. 10, we find that the left hand side is monotonically decreasing with E, and 
the right hand side linearly increases with E, therefore there can be at most one solution 
where l.h.s. = r.h.s., and thus only one bound state. In the case of Y(0) < 0, we also see 
that under the approximation of energy-independent coupling coefficients, the integrals over 
efe appearing in Eq. |TT| lead to terms of the form ln(ek)\ €k ={^ B T) ,E r Fermj y- Using this, along 
with recalling the coupling anisotropies A = T^/T^ and e = 1 — T^/Tp in order to define all 

: p — V-l — ojj. ,J- a — y\J- , J-p 



coupling coefficients in terms of T* = T: TS — (1 — e)T\ = AT; TP = (1 — e)XT, we find 



the following form for e, 



_ ^ i _ Fj^th ffb E 2, E 3 , E s Fermi , E d Fermi ) 



where F(Eq, E±, E 2 , E 3 , E Fermi , E Fermi ) is an algebraic function of all the energy parameters. 
The important feature to point out in Eq. [12] is the overall similarity with Eq. [8] The fact 
that a similar relationship between the coupling and anisotropy parameters of the system has 
been found from two very different approaches is remarkable: the creation of the relationship 
out of empirical analysis of the GME for the open system, and the analytical derivation from 
a bound state analysis of the eigenvalue equation. 

The formation of the bound state can be linked to an interference effect between trans- 
mission resonances due to the presence of several transport channels in the system — a 
feature eluded to in Fig. |4j Resonances of this type have been observed in diverse systems 
such as in quantum billiards,^ laser induced continuum structures in atoms,^ and quan- 
tum dotsP^ In these works, it has been revealed that the formation of bound states in the 
continuum is directly related to the shape of the system and anisotropy of the coupling 
parameters. This point is specifically reflected in our findings as we only observe the FSC 
resonance under coupling anisotropies. These coupling anisotropies may play the role of a 
path difference between the channels, where specific path differences may quench transport 
through the system in an effect akin to negative differential conductance.^ In such a case 
the FSC is enhanced — forming a resonance — since the dominant decoherence pathways 
are also effectively quenched. This may also help explain the presence of the single dominant 
resonance in the range of the orbital anisotropy in our model, since only a specific energy 
and coupling configuration will lead to transport quenching. 
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VII. CONCLUSIONS 



In conclusion, we have shown that the FSC of a fermionic system can be decoupled from 
the evolution of the population probabilities and the HSC of the system, even under the 
presence of boson-mediated relaxation and single-particle exchange with a bath. We demon- 
strated a resonant regime in this FSC where the decoherence times extend far beyond that 
of the HSC, and we have been able to relate these results to an analytically derived relation- 
ship in the coupling parameters of the system arising from the formation of a bound state 
between the system and reservoirs. We anticipate the primary difficulty in the experimental 
confirmation of our theory to be the establishment of the FSC in the first instance. What 
we have shown here is that such a FSC can exhibit a remarkable, robust, and long-lived 
resonance even in the presence of ordinarily very destructive perturbations such as particles 
tunnelling into and out of the system. A plausible configuration to obtain this FSC reso- 
nance effect is one where the particle reservoirs are macroscopic systems in thermodynamic 
equilibrium, as envisioned in the present work, and whose description is properly and most 
conveniently characterised by a particle density (not particle number) in a grand canoni- 
cal ensemble. Subsequent experimental and further theoretical efforts extending this work 
may shed light on the fundamental problem of decoherence, and on quantum information 
processing in a semiconductor environment. 
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